Surface pinning in amorphous ZrTiCuNiBe alloy 
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We have measured the amplitude and the phase of an electromagnetic (EM) field radiated from 
superconductor (amorphous ZrTiCuNiBe alloy) in the mixed state due to interaction of the flux 
lattice with an elastic wave. The results undoubtedly point to an essential contribution of a surface 
pinning into the flux lattice dynamics. We propose a model that describes radiation of EM field 
from superconductors with non-uniform pinning. The model allows to reconstruct the viscosity and 
the Labush parameters from the experimental data. The behavior of the Labush parameter can 
be qualitatively explained in terms of the collective pinning theory with the allowance of thermal 
fluctuations. 
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Soft type-II superconductors were studied intensively 
during last four decades with the goal to investigate var- 
ious aspects of vortex matter dynamics (see the compre- 
hensive review by Brandt pj). Nevertheless, many ques- 
tions that require further investigation remain. Among 
these questions is the problem of relation between the 
surface and the bulk pinning. Up to now the dynamics 
of the vortex state in such an exclusively non-uniform 
situation was studied by the surface impedance method 
2j. In this paper on the example of the amorphous 
Zr4i.2Tii3.sCui2.5NiioBe22.5 alloy we demonstrate the 
abilities of a new method based on excitation of the vor- 
tex lattice oscillations by a high-frequency sound wave. 

The essence of the method is the following. A super- 
conductor situated in lower half-space (z < 0) is sub- 
jected by a constant magnetic held H || z. A transverse 
elastic wave propagating along H and polarized in x di- 
rection produces transverse (with respect to H) oscilla- 
tions of the vortex lattice caused by pinning forces and 
viscous friction forces, and, consequently, induces electro- 
magnetic (EM) held. An antenna receives EM field (with 
E v and H x components) radiated through the elastically 
free surface of the sample (the surface perpendicular to 
the direction of propagation of the elastic wave) . While 
similar experimental setup was already applied for the 
study of type-II superconductors 0,0], the key new point 
is measuring both the amplitude and the phase of EM 
field (more accurately, the changes of these quantities). 

In the uniform case and in the local limit (q 3> I , 
q is the wave number and I is the mean free path) the 
components of EM field at z = are given by a simple 
expression 
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FIG. 1: Amplitude and phase of EM field measured for the 
sample with imperfect surface. Inset - behavior of \E\/H near 
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Hc2 (1 - experiment, 
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2 - computed from Eq. (7) with a 
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where u(z) — uq cos qze %u)t is the elastic displacement, k 2 
is the square of complex wave number of EM field in a 
conductor, and c is the light velocity. 



Eq. Q is applicable for normal as well as for super- 
conducting state of metal. In normal state k 2 = k\ = 
Airiujao/c 2 ((To is the static conductivity in normal state) 
is imaginary valued quantity. In superconducting state 
at small magnetic fields (H ~ H c \) k 2 = k 2 = \^ 2 (Aj, 
is the London penetration length) is real valued. In well- 
developed Shubnikov state k 2 = k 2 n = 4ir(iu>rj + oll) j H 2 , 
where 77 > is the viscosity parameter and > is the 
Labusch "spring" parameter. We count the phase ip of 
EM field from its value at small (H = H c \ + 0) magnetic 
field. As follows from Eq. (JTJ, the phase is positive at 
all H > H c i and it approaches <p n < 90° (the phase in 
normal state) at H — > H C 2- 

The initial goal of our experiment was to obtain exper- 
imentally the dependencies r)(H) and oil(H) using Eq. 
Q as it was done before for MgB2 .7] . We use the work- 
ing frequency uj/2it ~ 55 MGz. Details of the measuring 
procedure are described in ||. The results of our first 
experiments are presented in Fig. ^ 

Prior discussing these results let us give the values 
of parameters for the material under study that impor- 
tant for further analysis. The measured d.c. resistiv- 
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ity p ps 200 pfl-cm is practically temperature indepen- 
dent. The Alfer- Rubin effect (quadratic dependence of 
the sound attenuation coefficient on the magnetic field 
0) allows to determine the parameter = \q 2 /k 2 \ — 
91 ± 1, in excellent coincidence with the value of p 
measured (the sound velocity was found in 10]). The 
Hall constant, measured at room temperature, is posi- 
tive and very small: R H = ^ = (3.2 ± 0.4) • 10~ 25 
CGS units, that yields rather high density of the carri- 
ers n ps 2.2 ■ 10 23 cm -3 . Taking m = m e we obtain the 
relaxation time r ps 0.8 • 10~ 16 s and for the Fermi ve- 
locity vf ~ 10 8 cm/s we find I ~ 10 _8 cm, which is close 
to the dielectrization threshold. It was shown before [ll| 
that the alloy under study belongs to the family of weak 
coupling superconductors with the standard BSC energy 
gap A(0) fa 1.75T C (T c ps 0.85K). In the dirty limit the 
formula for the penetration depth |l2| can be rewritten 
as A^ 2 = fc 2 A tanh ^ , that yields Xl = 3 • 10~ 4 cm at 
T = 0.4K. The coherence length calculated from H C 2 is 
£(Q.4K)= 1.4- 10~ 6 cm and the Ginzburg-Landau param- 
eter is k = A L /£ ~ 200. 

As follows from Eq. , for (3 given above one could ex- 
pect (fin « 90° that contradicts with the result presented 
in Fig. [2(<^n ~ 120°). Let us show that the discrepancy 
found can be accounted for lowering of conductivity near 
the surface of the sample caused by imperfectness of the 
surface. 

We imply the following model dependence of the con- 
ductivity on z: a^(z) = ao v {l — p exp(z/z a )] (with p < 1 
and z a > 0). Then in normal state the electrodynamic 
equation for the EM field has the form (see |||): 



fE 



-^.-o(0S = o(C) cos(C), 



(2) 



where o(C) = a[l - p exp(C/Co)], a = k 2 nv /q 2 , k 2 nv = 
47riojao v I c 2 , (o = qz a . Here we use the dimensionless 
variables E — Ecj '{ilouqH) and C, = qz. The boundary 
conditions for Eq. (2) are that |2?(C)| is finite at all £ and 



C=o 



= 0. The latter condition is due to continuity 



dE/dC, 

of EM field on the conductor-vacuum interface and is 
valid with the accuracy S/Xem IS 10 -4 (6 is the skin 
depth and Xem is the wavelength of EM field in vacuum). 

The solution of Eq. J2J can be expressed through the 
Bessel functions of the 1-st kind of complex order v — 
2(oy/a on complex variable t(() = 2£o v / pa exp(£/2(o)- 
The field at the conductor-vacuum interface is deter- 
mined by the expression: 



E(z = 0) = 



dJ v [t{Q] 



C=o, 



dCJ«[t(QWQ cosC. 

(3) 

One can evaluate Eq. using the expansion of J v (t) 
in series in t and integrating each term analytically. For 
the case of interest it is enough to take into account first 
two terms of the expansion. 
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FIG. 2: Amplitude (solid curves 1,2,3) and phase (dotted 
curves la,2a,3a) of EM field measured for the sample with 
perfect surface at T = 0.4, 0.69, 0.8K, correspondingly (for 
T = 0.8K the results for the normal state are not shown). 
Circles - dependence \E(H)\ measured under decreasing of 
H . Inset - the same dependencies in the interval of small 
magnetic fields (all notations are the same as in the main 
plot). 



Analysis of Eq. © shows that the phase ip n increases 
under increasing p. In our case of extremely dirty con- 
ductor (\a\ ~ 10~ 2 ) and for Co ~ 1 we find that a rather 
small lowering of surface conductivity (~ 10%) compar- 
ing to the bulk one results in increasing of ip n up to 120°. 

Removing the layer ~ 50 /im from the surface of the 
sample by additional polishing with a fine powder (the 
size of grains ~ 1 2 /jm) we managed to eliminate the 
lowering of surface conductivity. The results of measure- 
ments (at various T) are shown in Fig. [3 One can 
see that in this case ip n ~ 90°. The main peculiarity 
of the dependencies observed is that at all temperatures 
the phase ip becomes negative at intermediate magnetic 
fields. This effect cannot be described by Eq. (1) under 
any reasonable variation of 77 and oll with H. 

In what follows we will argue that the effect observed 
can be explained by non-uniform pinning near the sur- 
face of the sample that takes place even in a situation 
with uniform conductivity (and, consequently, the uni- 
form parameter 77). Since peculiar behavior of the phase 
are observed in samples having less perfect as well as 
more perfect surface we consider the non-uniform pin- 
ning be an intrinsic property of the surface of the mate- 
rial under study. 

EM field in the mixed state is described by the system 
of equations, namely, the Maxwell equation, the matter 
equation, that determines the value of the current in the 
two-fluid model, and the equation of motion for the vor- 
tex lattice kj: 



(4) 



d 2 ~E Aniuj, 1 . 

-j x H + iw?y(u - u„) + a L {u - u v ) = 0, (5) 



where u,; is the displacement of the vortex lattice. Here 
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FIG. 3: Solid curves 1,2,3 - values of atL reconstructed from 
the data in Fig, with using Eq. 10; dashed curves la, 
2a, 3a - computed dependence a e /^(H) (see text) at T = 
0.4, 0.69, 0.8K, correspondingly. 



FIG. 4: Values of uirj reconstructed from the data in Fig, |5] 
with using Eq. J7J. Inset - the slop S of u)T) for intervals of 
linear increase of wrj with H (circles); solid line -linear fit. 



we neglect the normal component of the current. 

Non-uniform pinning can be modelled by z-dependent 
Labush parameter. We specify this dependence as 
= otLv + o>L S exp(z/z p ). For the magnetic field 
satisfied the inequality k 2 ^ \iojr) + cvl| we obtain 
from (JSJ) the equation for EM field that coincides in 
form with Eq. J2J) with the same boundary conditions 
and the quantity a(C) defined as 
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iuirj + a Lv 



■ exp ■ 



y 4tt y Air SU 



A 



Cexp-^, (6) 
Co 



where Co = Q z p- 

Taking into account first two terms in the expansion of 
the Bessel function we obtain from Eq. J2J the approxi- 
mate expression 
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Eq. JJJ can be used to reconstruct the dependencies r](H) 
and o.l{H) from the experimental data presented in Fig. 
13 if the parameters Co and aL V /&Ls are specified. These 
parameters are not known, but there are limitations on 
the choice of them. We require the functions r){H) and 
a ^(H) be smooth and real valued. Numerical analysis 
shows that these requirements satisfy for a^. u <C a^g 
(it means that the bulk pinning can be neglected) and 
Co < 0.5. 

The dependencies ctLs(H) obtained for Co = 0.5 are 
shown in Fig. |21 

There are physical reasons for Co to be close to Co = 0.5. 
One can see from Eq. 10 that smaller values of Co cor- 
respond to larger values of «l s . But down to the lowest 
temperatures available in our experiment we do not find 
(see Fig. |5J) any signs of freezing of the magnetic flux 
(the irreversibility line). It means that the pinning is 
quite weak and it is caused, most probably, by point de- 
fects. The estimates (see further) based on the collective 
pinning (CP) theory 01 show that the values of a^ s 



for Co = 0.5 (presented in Fig. |3J) are close to maximal 
possible value of this parameter. 

It is interesting to note that Co ~ 0.5 corresponds to 
z p ~ \l- Probably, it is a fingerprint of that the image 
forces are responsible, in some part, for a formation of 
z-profile of the pinning potential. 

The dependencies r)(H) are unsensitive to the choice 
of Co and linear in a rather wide interval of H with the 
slope S cx H C 2 (FigOJ. In this interval they have the 
form n(H) = (2 ± 0.05)HH c2 a n /c 2 . 

The CP theory 0] allows to compute the field and 
temperature dependence of the Labush constant using 
only one free parameter. This parameter is connected 
with the quenched disorder of the vortex lattice caused 
by a random pinning potential. It is convenient to choose 
a field H sv < 0.5H C 2 as such a parameter [li). This field 
separates the regime of single vortex pinning (SVP) from 
the regime of bundle vortex pinning (BVP) . In the SVP 
regime (H < H sv or H > H c2 — H sv ) the value of is 
given by the expression «l « CQQa~ 2 (H sv ), where C'qq = 

(8ttA1~P (l ~ ITa ) ^ s s ^ ear modulus of the vortex 
lattice, a v (H) = ^/Wq/H is the vortex lattice constant, 
and <f>o is the quantum of magnetic flux. In the BVP 
regime (H sv < H < H c2 — H sv ) the Labush parameter 
is ai w CqqR~ 2 (H), where the collective pinning radius 
R c is given by the formula |l5l | 



R c (H)^a v (H)exp- 



H(H c2 - H) 
H sv (H c2 — H sv ) 



- 1 



(8) 

The expressions presented determine the dependence 
C(l(H) almost symmetric relative to the H = 0.5H c2 
line with the shape varied from the bell shape at H sv ~ 
0.5H c2 to the double- hump shape at H sv <C H C 2- 

One can see from Fig. [3]that while the value of ob- 
tained from the experimental data is of the same order 
as one given by the CP theory estimates, but the de- 
pendence aL(H) is strongly asymmetric relative to the 
H = Q.5H C 2 line. We consider that such behavior of 
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are caused by thermal fluctuations. 

To evaluate the effect of thermal fluctuations one 
should add a random force fa into Eq. © 1(|. Strictly 
speaking, the returning force jH/c = (H 2 / '4ir)d 2 u v / 'dz 2 
in Eq. (JSJ has to be found from joint solution of the sys- 
tem (0J, but for the estimates one can replace d 2 /dz 2 
with — q 2 . In this case Eq. (0 can be rewritten in the 
form of equation of diffusion of the Brownian particle 

^ = -^(V° + Vu + V p )+fa, (9) 

where w = u — u v , the "tildes" indicate that the cor- 
responding quantities are given per one diffusing "parti- 
cle": y = y^QL c /H (L c is the collective pinning length), 
V = {q 2 H 2 /£,Tr)w 2 , V u = {q 2 H 2 /4ir)uw, and V p is the 
pinning potential. The latter is modelled by a three-well 
potential 

r ( w + d) 2 w <~d/2 

V P = ^-lw 2 \w\ < d/2 (10) 

1 [ (w - d) 2 w> d/2 

where d is the distance between minima of the pinning 
potential. The linearized Fokker-Planck equation that 
corresponds to Eq. © with the pinning potential lfTT)|l 
has the exact solution in terms of hypergeometric func- 
tions. At ai/(q 2 -^r) <C 1 the following simple esti- 
mate for the averaged displacement of the vortex lat- 
tice is found (u v ) ~ u ■ An^iujri + a e / ^ ) / (q 2 B 2 ) , where 

ag* = a L (l - ^ce" c2 ) and c 2 = ''y d 2 . One can 
see from the comparison of this formula with Eq. 
that the thermal fluctuations may result in an essential 
reduction of the effective Labush parameter remaining 
the viscosity parameter unchanged. For general case the 
dependencies a e / (H) obtained from the solution of the 
Fokker-Planck equation are shown in Fig. [3] To achieve 
semi-quantitative agreement between the theoretical re- 
sults and the experimental data at T — 0.4K we choose 
H sv as 0.15H C 2 and d « 5 • 10~ 7 cm. Other important 
points of the fitting procedure are the following: I)To 
agree the maximum value of ot^ (at 0.4K) with the ex- 
perimental one we take for that enter into equation 
for C 66 the value X L (0AK) = 1.2 • 10" 4 cm. The stan- 
dard temperature dependence of A^ 0| is implied. II) 
We set L c — a v (H) . Ill) The quantity d is assumed tem- 
perature dependent (d(T) cx IV) We imply that 
temperature dependence of H sv is determined by the 8T C 
pinning: H sv /H c2 oc (1 - (T/Ty 2 )" 1 / 3 0. Such an as- 
sumption looks quite reasonable if one takes into account 
that there is a superconducting phase with higher T c on 
the surface of our sample [llj] . 

One can see from Fig. [21 that computed temperature 
and field dependencies of az, are in qualitative agreement 
with the experimental results. 

It is interesting to note that if we substitute into 
Eq. (0) we obtain \E\/H that has a local minimum near 



H C 2- Such a minimum is observed sometimes in our ex- 
periment (Fig. 2]) as well as in Q. Physically this min- 
imum is connected with the screening of EM field radi- 
ated from dipper regions of the sample by the surface 
layer with reduced penetration depth. Reminiscence of 
this minimum is also present in Fig. [21 under transition 
from the normal to mixed state the increase of \E\ begins 
only after substantial lowering of if. 

In conclusion, we propose the new method of investi- 
gation of dynamical characteristics of the vortex matter 
that consists in measuring the amplitude and the phase 
of EM field radiated from the sample under excitation 
of the vortex lattice oscillations by elastic wave. It is 
established that unusual behavior of the amplitude and 
the phase of EM field are accounted for the surface pin- 
ning. The parameters measured (viscosity coefficient and 
Labush parameter) are in good agreement with the es- 
timates obtained from the theories of vortex matter dy- 
namics. 
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